The electrical double layer for a fully asymmetric electrolyte around a spherical 

colloid: an integral equation study 



O 

o 

(N 
o 

Q 



o 



o 



G. Ivan Guerrero-Garci'a,^ Enrique Gonzalez- Tovar,^''^ Marcelo Lozada-Cassou,^ and F. de J. Guevara-Rodriguez^ 

^ Instituto de Fisica, Universidad Autonoma de San Luis Potosi, 
Alvaro Obregon 64, 78000 San Luis Potosi, S.L.P., Mexico 
^ Programa de Ingenieria Molecular, Instituto Mexicano del Petroleo, 
Eje Central Ldzaro Cardenas 152, 07730 Mexico, D.F., Mexico 
(Dated: February 2, 2008) 

The hypernetted chain/mean spherical approximation (HNC/MSA) integral equation is obtained 
and solved numerically for a totally asymmetric primitive model electrolyte around a spherical 
macroparticle. The ensuing radial distribution functions show a very good agreement when com- 
pared to our Monte Carlo and molecular dynamics simulations for spherical geometry and with 
respect to previous anisotropic reference HNC calculations in the planar limit. We report an analy- 
sis of the potential vs charge relationship, radial distribution functions, mean electrostatic potential 
and cumulative reduced charge for representative cases of 1:1 and 2:2 salts with a size asymmetry 
ratio of 2. Our results are collated with those of the Modified Gouy-Chapman (MGC) and unequal 
radius Modified Gouy-Chapman (URMGC) theories and with those of HNC/MSA in the restricted 
primitive model (RPM) to assess the importance of size asymmetry effects. One of the most striking 
characteristics found is that, contrary to the general belief, away from the point of zero charge the 
properties of an asymmetric electrical double layer (EDL) are not those corresponding to a symmet- 
ric electrolyte with the size and charge of the counterion, i.e. countenons do not always dominate. 
This behavior suggests the existence of a new phenomenology in the EDL that genuinely belongs to 
a more realistic size-asymmetric model where steric correlations are taken into account consistently. 
Such novel features can not be described by traditional mean field theories like MGC, URMGC or 
even by enhanced formalisms, like HNC/MSA, if they are based on the RPM. 

PACS numbers: 61.20.-p, 61.20.Gy, 61.20.Ja, 61.20.Qg. 
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I. INTRODUCTION 

The electrical double layer (EDL) is the structure 
formed by electrolyte ions around a charged surface, usu- 
ally that of a colloid or electrode. An understanding of 
the EDL properties is a crucial matter for science and 
technology because of the large variety of related ap- 
plications, that range from colloidal stability, electroki- 
netics and the description of biological systems to daily 
manufactured products as inks, paint emulsions, foods or 
medicaments • As an illustration, the EDL deter- 
mines completely the value of the zeta potential of a col- 
loid in electrophoretic motion 0, Q . The zeta potential, 
which is directly related to the measured mobility in elec- 
trophoresis experiments, is a central quantity in colloid 
science and many standard techniques of characteriza- 
tion, stabilization and separation of colloidal suspensions 
rely critically on its knowledge 5, 6, 7, 8]. Depending on 
the form of the dispersed macroparticle (or electrode) 
is that we are dealing with an EDL of particular geom- 
etry, e.g. planar, spherical or cylindrical. The planar 
case is, by far, the most studied instance, however, the 
spherical EDL (SEDL) deserves special attention for its 
obvious relation to dispersions of globular proteins, mi- 
celles, polymer beads, dendrimers or many other nearly 
spherical organic or inorganic macroions Qi&^^J}- 

Since the emergence of this topic, the unquestionable 
relevance of the EDL has been paralleled by an intense 
search of an adequate theoretical description of such 



charge distribution. In order to build a successful theory 
of the EDL, a good model is essential. Among some more 
refined representations, nowadays, the usual election for 
the EDL is the minimal but able restricted primitive 
model (RPM) . In this well-known idealization of a single 
EDL a "colloid" or "electrode" is mimicked by a rigid and 
uniformly charged object, which is immersed in an elec- 
trolyte solution constituted by a structureless solvent me- 
dia and various ionic species, where the latter are treated 
as equally-sized hard spheres with punctual charges at 
their centers. To avoid image effects it is customary to 
suppose a uniform dielectric constant permeating all the 
space. During the past decades the RPM-EDL has been 
studied comprehensively via modern statistical mechan- 
ics approaches, namely integral e quat ions 0, 0, 
111 111 HI El Ip, 1?, 'HEljlpJIi , density func- 
tional theories 0, |25., 26, 2j|2a [Mm 32 33| and 
mean electrostatic potential schemes [SJ, IH, IH • Ad- 
ditionally, the usage of simulation techniques to evaluate 
the properties of these coulombic systems has supplied a 
large body of "experimental" data useful to test the exist- 
ing theories [13, IMIM 113 El El El El El El El- At 

the end, and on spite of its apparent simplicity, all the 
collected evidence sustains the RPM as a concise and 
thriving model which, by embodying the most impor- 
tant interactions in EDL systems, i.e. electrostatic and 
hard-core potentials, is capable to provide not only an 
essentially correct account of the associated thermody- 
namic properties but also to predict characteristic phe- 
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nomena of the EDL such as layering, charge inversion, 
charge reversal and overcharging. Briefly, we remember 
to the reader that the stratified (oscillating) ionic dis- 
tributions next to a colloid arc referred to as layering. 
On the other hand, charge reversal means the excessive 
compensation of the native colloidal charge, prompted 
by strongly attracted counterions, leading to an effective 
macroparticle's charge of reversed sign. This counteri- 
ons layer, in turn, promotes a layer of colons adsorption 
which, consequently, produces charge inversion, i.e. a 
layer of inverted charge. Complementarily, overcharg- 
ing is the unusual adsorption of colons to the surface of 
the colloid, which increases its original charge In 
fact, the continued progress of theories founded on the 
RPM has made of them a tool of quantitative capacity 
to interpret experimental results. As an example, very 
recently, the new Primitive Model Electrophoresis ap- 
proach '4^ [H^l has incorporated consistently the ionic 
size effects into the electrokinetic equations and forged 
an enhanced treatment of electrophoresis, which fits cor- 
rectly measurements in multivalent electrolytic ambients 
and explains reversed mobilities [43, 113, IsH ■ 

Historically, the RPM-based theoretical attempts 
surged as to better the classical Poisson-Boltzmann (PB) 
portrait of the EDL. In the conventional PB treatment, 
instead of hard-sphere ions, a punctual electrolyte is as- 
sumed and the ionic distributions are obtained by solv- 
ing the Poisson-Boltzmann differential equation, thus 
eliminating the very important interionic steric corre- 
lations present in the RPM. Although the assumption 
of point-ions could be acceptable in some cases (e.g., 
if low charged colloids in very diluted univalent salts 
are being considered), for high-coupled systems the in- 
adequacy of such h ypq tesis has been clearly attested 
[1^ Il3l IT^ l3^ I53L1m|. As a consequence, the orig- 
inal bare PB theory is unable to predict any of the 
nonmonotonic characteristics mentioned in the previous 
paragraph. Since its early stages, the patent faults of 
the PB equation were already noticed [sM Is^ Is^ Isl ] 
and, accordingly, a first modification to the simple PB 
treatment was devised by introducing the concept of the 
Stern layer or the Helmholtz plane. In this supplemented 
model, the only inclusion of a unique distance of clos- 
est approach between the ions and the macroparticle re- 
sulted in the normalization of the otherwise unphysical 
values of the radial distribution functions (RDFs) close 
to the contact and, eventually, extended the applicabil- 
ity of the PB-Stern or modified Gouy-Chapman theory 
(PB-S or MGC) to 1:1 systems with moderate charge or 
concentration j^ls^- Nevertheless, MGC remained ex- 
empt of nonmonotonical features and, thus, unsuccessful 
to cope with more demanding conditions (viz. multiva- 
lent and/or highly concentrated electrolytes or, else, high 
surface charges). 

Up to now, most of the theoretical work on EDL has 
supposed the valence as the unique source of ionic asym- 
metry. Obviously, in nature this is not the situation and, 
out of the many possible sophistications to the model 



(e.g. a discrete solvent, more species, etc), the con- 
sideration of ions with different sizes is a first choice. 
Some works have been already published along these lines 
[63,|6ll|6l|6l|6l|6lE^,'6^ but, notably, the pioneer- 
ing attempt to include distinct ionic sizes was formu- 
lated inside the framework of the PB equation by Val- 
leau and Torrie who used Stern layers of unequal 
extent for the counter- and colons. It must be noted 
that in this unequal radius MGC (URMGC) scheme the 
ions are, at the same time, voluminous and punctual 
objects, i.e., an electrolytic particle behaves as a hard 
sphere when interacting with the colloid (being dissimi- 
lar the colon and counterion distances of closest approach 
to the macroparticle), but among ions they "see" each 
other as points. In summary, such addition to MGC was 
rewarding since exposed interesting effects not observed 
in the size-symmetric case, e.g. the displacement of the 
potential of zero charge (PZC), the occurrence of oscilla- 
tions in the radial distribution functions and the apparent 
dominance of counterions to determine the properties of 
asymmetric EDL systems. At first, all these phenomena 
seemed appealing since, for instance, the non-zero PZC 
allowed the possibility to interpret experimental data by 
obviating specific adsorption, and the occurrence of fluc- 
tuating URMGC ionic densities and electrostatic poten- 
tials meant that a "point-ion" theory should not be nec- 
essarily uniform. Furthermore, for a n:n electrolyte with 
a common counterion size, the convergence between the 
MGC and URMGC outcomes was so clear that founded 
the now widely accepted fact that the counterion size is 
what really matters for the EDL. Quoting Valleau and 
Torrie: "...we expect the double layer properties of a di- 
lute (asymmetric) electrolyte to become similar to those 
of a completely symmetric electrolyte having an effective 
size equal to that of the counterion. (This remark will 
be asymptotically exact for large fields in the Poisson- 
Boltzmann theory)..." |60j |. Notwithstanding, due to the 
absence of steric interionic correlations, all the new fea- 
tures in URMGC are constrained in space and can be 
rationalized purely in terms of energetic arguments, leav- 
ing out determinant entropic contributions 48, 68. .69i| , so 
that, at the end, URMGC will be incapable to produce a 
full description of the phenomenology happening in gen- 
eral size-asymmetric systems. Even more, we will make 
evident that an analysis based on the PB viewpoint of a 
totally asymetric EDL not only is partial but, regarding 
the counterion predominance, can lead in wrong direc- 
tions. 

For all the stated, to improve the theory of EDL the 
asymmetry in ionic size must be included in a more con- 
sistent way. That can be accomplished by means of the 
unrestricted primitive model (PM) of an electrolyte, for 
which the condition of equal size and charge in the RPM 
is relaxed. Not surprisingly, in the last years, the PM- 
EDL for plane electrodes has been studied through the 
avenues of integral equations and the mean electrostatic 
potential 64, 65, 66]. Amid them, the Greberg and Kjel- 
lander's paper 66] represents a valuable contribution for 
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its detailed appraisal of the effects of size asymmetry 
on the ionic structure and, chiefly, for its thorough dis- 
crimination of the different contributions to charge in- 
version and to the mean force. In any case, the body 
of existent literature on PM-EDL shares two limitations, 
namely they have dealt exclusively with the planar ge- 
ometry, ignoring the undoubtfully notable spherical in- 
stance JO], and have focused on the structural proper- 
ties or limited themselves to describe semiquantitatively 
the behavior of the potential versus charge relationship. 
Besides, and most importantly, due to the conditions ex- 
plored, the majority of modern studies going beyond PB 
[Ulli m HI m m IzJ have subscribed the wonted 
conclusion that the counterion-macroparticle interaction 
is what determines the EDL properties: ^^...Away from 
the pzc, the double layer is dominated by the counteri- 
ons and, for a binary system at least, the properties of 
the double layer are essentially those of a symmetric salt 
with the charge and diameter of the counterion...^^ |63l |. 
We will show that this is not always the case. 

Consequently, in the present paper we report a com- 
prehensive investigation of the PM-EDL for a spheri- 
cal model colloid using integral equations and simulation 
techniques with a two-folded aim: i) to produce a descrip- 
tion of the significant PM-SEDL system using a trustable 
statistical mechanics treatment and ii) to reivindicate the 
importance of size asymmetry in EDL studies in order 
to get a more faithful picture of a model colloidal sus- 
pension and as the origin of an enriched phenomenology, 
previously unseen in the symmetrical conditions. The 
rest of the paper is organized as follows. In Section II 
(and Appendices) we describe the model and give the 
details of the HNC/MSA theory for the PM-SEDL and 
of the employed simulation methods. The results, com- 
parisons and their discussion are included in Section III, 
and we close in Section IV with our conclusions and some 
prospectives for future work. 



II. THEORY AND METHODS 

A. The Basic Model 

Our representation of the SEDL is constituted by a 
rigid, charged spherical colloid of diameter D and surface 
charge density ctq i surrounded by a continuum solvent of 
dielectric constant e. The macroion is in contact with 
two ionic species which are treated as hard spheres of 
diameter Ri {i — 1,2) with embedded point charges of 
valence Zi at their centers. It is assumed for definiteness 
that R2 > Ri- The interaction potential between the 
macroion and an ion of type i is then given by 



00, 



r < 



D+Ri 



with e the protonic charge. In turn, the interionic poten- 
tial is 



00, 



r < 
r > 



Ri+Rj 
2 ■ 

Ri+Rj 



(2) 



The Stern layer or, more properly, the Helmholtz 
plane, already mentioned in the introduction, is the geo- 
metrical place corresponding to the closest approach dis- 
tance for the electrolyte ions to the colloid. If we consider 
an electrolyte formed by a pair of ionic species of unequal 
sizes, the closest approach distance to the surface for the 
smallest component, {D + i?i)/2, determines the inner 
Helmholtz plane (IHP), whereas the corresponding dis- 
tance for the largest species, {D + i?2)/2, establishes the 
outer Helmholtz plane (OHP). In the limit of identical 
sizes the IHP and OHP coincide and the usual definition 
of the Helmoltz plane is recovered. 



B. The HNC/MSA equation for the PM-SEDL 

In general, the Ornstein-Zernike equation and the hy- 
pernetted chain closure (HNC) for a multicomponent 
mixture of M species are 

M „ 

hij {r 12) ^ Cij{r 12) +'^pi hu{ris)cij{r32)dV (3) 



and 



Cijir 12) = -l3Uij(r 12) + hij{r 12) - ln{hij{r 12) + 1), (4) 

i,j = l,2...M, 

where hij{r 12) is the total correlation function for parti- 
cles 1 and 2 of species i and j, distant by ri2, which is 
related to the radial distribution function by gij{ri2) = 
hij{r 12) + 1, Cij(?'32) is the direct correlation function for 
the pair 2 and 3, pi is the bulk number concentration 
of each one of the species, C/y (ri2) is the direct interac- 
tion potential, dV = d^ria is the differential volume, and 
(3 = l/{kBT), where fc^ is the Boltzmann's constant and 
T the absolute temperature. 

Applying the Direct Method 7?, "t^, Ti^ , let the species 
M correspond to macroions at infinite dilution in a binary 
electrolyte. Then Eqs. ^ and for species M and j 
become 



hMj{ri2) = CMj{ri2) + '^Pi J hMiiri3)cij{r32)dV (5) 



1=1 



and 



(1) CMj{ri2) = -PUMj{ri2) + hMj{ri2) - ln{hM]{ri2) + 1), 

(6) 
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i = l,2. 

If the HNC closure, Eq. ©, is used for the macroion- 
ion direct correlation function in the Ornstein-Zernike 
equation |(SJl, we get 



0, 



r < 



9Mj{r) 



exp{~(3UMj{r) + flMjir)} , r > 2^., 



with 



1=1 •' 



(7) 



(8) 



where r = ri2, t = ri3 and s ~ r23 =| t — f \ are the 
distances between particles 1 and 2, 1 and 3, and 2 and 
3, respectively, and dV = senO d9 dcj) dt is the volume 
element in spherical coordinates. The geometry of the 
system is illustrated in Fig. ^ 

It must be noted that the sizes of the macroion and of 
the ionic species are arbitrary, and the only demanded 
condition is the electrolyte bulk electroneutrality: 



Pizi + P2Z2 = 0. 



(9) 



In respect to the ion-ion direct correlation functions, 
Cij(s), in the bulk MSA theory they can be written as 



Cy (s) 



(10) 



These functions are dependent on the size and charge of 
the ionic species and on other parameters of the system 
such as the temperature and the dielectric constant, and 
are stated explicitly in the Appendix A. 

Introducing the MSA direct correlation functions in 
Eq. we arrive to the following nonlinear system of 
integral equations (the index M is dropped because there 
is only one macroion present): 



5i(r) = exp{-l3Ui{r)+ pi / gi{t)cii{s)dV 



Pi / Cii{s)dV + p2 / g2{t)c2i{s)dV - 
D + Ri 



P2 j C2l{s)dV 



r > 



(11) 



and 



g2ir) = exp\-pU2ir)+pi / giit)ci2is)dV 



Pi J Ci2{s)dV + P2 J g2{t)c22{s)dV- 

P2 C22(s)dV}, r> . 



(12) 



From the specific Cij (s) definition and considering the 
electroneutrality restriction, Eqs. (|11|) and H12|l can be 
recast in a more compact form as 



gi{r) = exp {lo^r) + pili(r) - pihir) 
+P2h{r) - P2h{r)} 



and 



92{r) 



with 



and 



^zie47r(f )^cro fl 
ksTe 



Ii{r) = / gi{t)cii{s)dV, 



Hr) = / cii""(s)dy. 



h{r) = / g2{t)c2i{s)dV, 



h{r) = / C2i™*(s)dy, 



Mr) 



2:2e47r(f ) ctq fl 
ksTe 



Hi{r) = / giit)c,2{s)dV, 



Hiir) = / C22™'(s)dK 



(13) 



exp{Ho{r) + piHi{r) - piH2{r) 
+P2H3ir) ^ P2Hi{r)} , (14) 



(15) 



(16) 



(17) 



(18) 



(19) 



(20) 



(21) 



H2{r) = I ci2'"*{s)dV, (22) 



H:,{r)^ J g2{t)c22is)dV, (23) 



(24) 



The detailed form of these terms is given in the Appendix 
B. 
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As a particular case, when Ri = R2 in Eqs. 
and ©I the equations l(T^ and ffT^ are reduced to the 
HNC/MSA expressions for the RPM-SEDL, deduced by 
Gonzalez- Tovar and Lozada-Cassou 16]. If the interionic 
size correlations are neglected in the diffuse EDL but the 
closest approach distances between the macroion and the 
ionic species are conserved, by setting i?i = i?2 = in 
Eq. |(2Jl and letting i?i and R2 arbitrary in Eq. it 
is found that h^r) = h{r) = H2ir) = Hi{r) = and 
the integral version of the URMGC theory is obtained. 
Additionally, if i?i = i?2 = in Eq. Q and i?i = R2 
in Eq. the HNC/MSA equations are reduced to the 
corresponding MGC equations in integral form. 

Derived from the ionic profiles, two quantities of spe- 
cial interest are the mean electrostatic potential (MEP) 
and the integrated surface charge density, which, consid- 
ering the electroneutrality condition, can be written as: 



, , , Aire f°°y-^ , ^ 

^ J r 



t 

r 



dt 



(25) 



and 



cr(r) 



e 



/oo r 



fdt, 



(26) 



respectively. Such quantities yield important parameters 
of the system when evaluated at specific positions. For 
example, if r = D/2 Eqs. (|25|l and (|26|l become the 
MEP at the macroion's surface, ipo, and the native sur- 
face charge density of the macroion, cfq, respectively: 



associated with the MEP at the closest approach dis- 
tance between an electrolyte and a charged surface 
in EE III III llllli- Following this convention, in a 
size- asymmetric or PM electrolyte ^ — ^ihPi and when 
the ions are equally-sized (in the RPM) ^ = ipiHP = 

TpOHP- 

The integrated surface charge density is a very conve- 
nient magnitude for infinite surfaces such as unbounded 
planes or cylinders. Nevertheless, for finite surfaces, as 
spheres, ellipsoids or spherocylinders, is useful to consider 
instead the cumulative reduced charge (ORG): 




t^dt , (31) 



where Q{r) = ATTr'^a{r). Beyond the macroparticle's sur- 
face this quantity corresponds to the total charge (native 
plus adsorbed) inside a sphere of radius r, normalized 
with the bare charge over the colloid. If D/2 < r < 
{D + Ri)/2 the GRC is equal to one, whereas for r 00 
this quantity goes to zero because of the electroneutrality 
condition. The GRC has also the property of indicating 
charge reversal when Q*{r)Qo < 0. Moreover, charge 
reversal and overcharging are conveniently quantified by 
Q*(r) as a fraction or multiple of the native bare charge 
compensated locally up to certain distance. 



C. Computational methods 



■00 



Aire 
e 



+2 1 



dt 



(27) 



and 



rdt. 



(28) 



Furthermore, ior r = {D + Ri)/2 Eq. H25() corresponds to 
the MEP at the IHP, ^/jihp, whereas for r = (£> + i?2)/2 
the MEP at the OHP, i/jqhp, is obtained, i.e. 



D+Ri 



dt (29) 



and 



47re /"^ V- . N 

i'OHP^— j^}_^9^{t)p^ 



t' 



D+R2 



dt. (30) 



A fundamental quantity in electrokinetics is the zeta 
potential, (. This potential has been conventionally 



For a given value of ctq, the set of nonlinear integral 
Eqs. H13|) and (|14|l was numerically solved via the suc- 
cessive substitution (or Picard) method. In the past, 
this technique has been successfully applied to work out 
integral equations associated to similar coulombic fluids 
problems (12, j32j . We have produced accurate results 
by continuing the iterative procedure until the Euclidean 
norm between consecutive input and ouput of the RDFs 
was less than some prescribed small number, guarantee- 
ing an error inferior to 1%. As an extra test of con- 
sistency, at the end of the Picard cycle, the final ionic 
profiles were used to check to what extent they fuUfiUcd 
the total electroneutrality condition in the form of Eq. 
(|28|l . In all our calculations o-q was recovered up to a 
99.9 % at least. 

For comparison intentions, the static properties of 
some of our PM-SEDL samples were also computed by 
performing Monte Carlo (MG) and molecular dynamics 
(MD) simulations. In particular, we focused on 2:2, 0.5 M 
electrolytes with size asymmetry ratio of 2 (characterized 
by the parameters listed in Tabled bathing a macroion. 
This macroion was considered as a hard sphere of diam- 
eter D = 10 A with a charge ZmS located at its center, 
with Zm= 8 and e the protonic charge, which correspond 
to a surface charge density ctq = 0.407 C/m^. The solvent 
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enters the model as a uniform dielectric with an aqueous 
value of e = 78.5 at a temperature T — 298 K. In order 
to avoid image potential effects the dielectric constant of 
the macroion and the ionic species was considered equal 
to the solvent. The simulations were carried out within 
a cubic box with one macroion fixed in the center, sur- 
rounded by the electrolyte. The number of ions for each 
species was adjusted to satisfy the electroneutrality con- 
dition over all system: NiZie + N2Z2e + ZMe=0- The 
usual periodic boundary conditions with the minimum 
image convention for the ions were imposed to the simu- 
lation box, whose length (L) was considered sufficiently 
large to neglect macroion-macroion interactions and to 
obtain stable profiles (the specific parameters used are 
listed in Table P). Following Degreve et al. [4(1. l4ll|. 
the use of such a sufficiently large simulation box for 
SEDL systems owning a plasma parameter [t^IsoIIsiLIs^ 
Tp < 10 makes unnecessary a treatment a la Ewald of the 
coulombic interactions. Precisely, in our study Tp = 3.1. 

The MC simulations were performed in the canonical 
ensamble using the standard Metropolis method |83.l83|. 
The macroion-ion and ion-ion interaction potentials are 
given by Eqs. ^ and ||2Jl, respectively. In the thermal- 
ization process were carried out 4 x 10^ MC cycles to get 
the system into equilibrium and 2 x 10^ MC cycles were 
performed in order to take the canonical average after 
the equilibration. The aceptance ratio used was 0.4. 

In the MD simulations the numerical integration of the 
motion equations was accomplished through a reversible 
MD technique ^85] . This method ensures numerical sta- 
bility and is equivalent to the velocity Verlet algorithm. 
In order to produce data in a feasible way and to avoid 
the problem of an indefinable force at contact, a very 
steep but continuous interparticle potential in lieu of a 
hard core one was employed, i.e. 



A 



0, 



r > Ri. 



(32) 

where A = ksT, I = 30 and Rij is the closest approach 
distance between any pair of particles either macroion- 
ion or ion-ion. In consequence, the total potential among 
charges is the sum of Wij(r) and the Coulomb contribu- 
tion as defined in Eqs. ^ and The time step used 
was At — 10^'^r, where r — y^moRf/ {ksT) is the unit 
of time and mo is a characteristic unit of mass. The iner- 
tial masses of the components of the system are mi = mp, 
7712 = 8.0 77io and mmacroion — 13.03 777o. Finally, the 
formal runs involved more than 8 x 10*^ MD steps after 
equilibration. 



III. RESULTS AND DISCUSSION 

In order to assess the reliability of the HNC/MSA 
equation for the PM-SEDL a comparison of our re- 



sults with simulation data and, possibly, against alter- 
native theoretical approaches is called for. In the first 
case, collating with computer experiments is very use- 
ful to discriminate the adequacy of an approximate the- 
ory for a given model system. Additionally, confronting 
HNC/MSA with another available and, supposedly, more 
robust formalisms allows us to test its degree of accuracy 
and economy. 

Previously, it has been shown that HNC/MSA agrees 
very well with computer data for the RPM of pla- 
nar, spherical and cylindrical EDLs with symmetrical or 
asymmerical valences [HIIlliSlillliSliHlilii. On 
these grounds, now, we have opted for a certainly strict 
comparison, namely for a divalent electrolyte with ap- 
preciable size asymmetry. The two available options for 
dissymmetry were considered: counterions smaller than 
colons and viceversa. 

At this point, it is important to mention that even if 
the equations developed in the last section are sufficiently 
general to include asymmetry in size and charge, here we 
will restrict our attention to n:n electrolytes (i.e. | z+j = 
\ Z-\) to centre on the effect of different sizes. Following 
this convention, thereinafter, when we refer to PM we 
will only imply distinct ionic diameters. 

In Fig. 13 (main panel) we present the radial distri- 
bution functions for HNC/MSA and MC simulations for 
a macroion of diameter D = 10 A and surface charge 
density ctq = 0.407 C/m^ in a 2:2, 0.5 M electrolyte of 
diameters i?_ = 4.25 A and R+ = 8.5 A. The inset of 
Fig. 121 shows RDFs for the same system but collating 
HNC/MSA with MD data. The Fig. El contains a similar 
comparison between MC and MD simulations and the 
HNC/MSA theory for the same system considered in the 
Fig. 121 but now for a salt with interchanged ionic sizes, 
R- = 8.5 A and i?+ = 4.25 A. From these evidences, it is 
rewarding the excellent concordance between theory and 
simulation for the two possible asymmetries. 

Complementarily, we also made some comparisons be- 
tween the spherical HNC/MSA in the limit D ^ oo and 
published data of the anisotropic reference HNC theory 
(ARHNC) for a charged planar wall '66"|. We resorted to 
this limit because, to the best of our knowledge, there are 
no alternate structural results for a charged macrosphere 
in a size-asymmetrical electrolyte 'g^. Furthermore, at 
least for the planar geometry, ARHNC is an accurate de- 
scription that considers the inhomogeneous correlations 
in the electrolyte due to the wall and that fits very closely 
the simulations of the RPM-EDL. Thus, in Fig. ^(main 
panel) are displayed the HNC/MSA ionic profiles for a 
huge macroion of diameter D = 1000 x i?_ and charge 
density ctq — 0.267 C/m^, surrounded by a 1:1, 1 M elec- 
trolyte with i?_ = 4.25 A and i?+ = 8.5 A, and the 
ARHNC RDFs corresponding to a charged wall bearing 
the same ctq and under identical electrolytic parameters. 
To cover a more exigent situation, in the inset of Fig. 01 
we plot the HNC/MSA distributions of a 2:2, 1 M elec- 
trolyte of ionic dimensions i?_ — 4.25 A and R+ = 6.375 
A around a large colloid, characterized by D = 1000 x i?_ 
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and (To = 0.267 C/m^, along with the ARHNC out- 
comes for the same electrolyte in contact with a wall 
of equal ctq. Clearly, for both the uni- and divalent sys- 
tems, the coincidence between HNC/MSA and ARHNC 
is notable, considering that HNC/MSA is a simpler and 
less computationally demanding theory. Summarizing, 
the comparisons with Monte Carlo, molecular dynamics 
and ARHNC prove the trustabihty of HNC/MSA for the 
PM-SEDL and gives us the confidence to explore its pre- 
dictions in a wider range of conditions. 

A global and concise manner to analize the properties 
of the EDL for a large set of different states is in terms 
of the relationship between the MEP at some some point 
and the surface charge. In this respect, the first charge 
derivative of the MEP at the surface, is relevant 

because of its connection with the differential capacity 
m m El m S m m, whereas the MEP at the 
IHP, ipiHP, is conventionally associated with the zeta 
potential, a key quantity in electrokinetic phenomena. 
Besides, and as it will be shown, it is precisely in the 
behavior of these MEPs that the importance and exten- 
sion of the size asymmetry effects under study are evi- 
denced. Because of this, we have calculated the functions 
■00 (co) and V'/-f/p(co) for a size- asymmetric SEDL using 
the following four theoretical schemes: URMGC, MGC, 
HNC/MSA in the PM (HNC/MSApAf), and HNC/MSA 
in the RPM (HNC/MSA^pm). In a series of prelimi- 
nary runnings we have noticed that the manifestations 
of size asymmetry are more emphatic when the counte- 
rions are smaller than colons, thus, below we shall con- 
centrate on the case where the diameter of the coions 
is twice that of the counterions. For comparison pur- 
poses, the parameters of the systems are adjusted such 
as for URMGC and HNC/MSApm the unequal coUoid- 
counterion and coUoid-coion closest approach distances 
are shared, and for MGC and HNC /MS A/jpm the unique 
closest approach distance corresponds to that employed 
in HNC/MSApM for colloid- counterion. Therefore, in 
the rest of our calculations when the ionic diameters of a 
PM electrolyte, i?_ and are stated it is meant that 
the different theories are solved for the following elec- 
tions of the colloid-ion and ion-ion distances of closest 
approach, di and dij, respectively (with ctq > 0): 



and di 



= d.= 



D+R- 



forHNC/MSApM 
and URMGC, 

for HNC/MSAppM 
and MGC, 

(33) 

, for HNC/MSApM, 
for HNC/MSAppM, 
for URMGC and MGC. 

(34) 

For the theories mentioned in the foregoing the 
dependence of ■00 and ^ihp on the surface charge 
(To is given in Figs. Ela) and|SIb), respectively, for a 
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macroion with D = 160 A and (To > and a 1:1, 1 
M electrolyte with i?_ = 4.25 A and R+ = 8.5 A. In 
those figures, for any non-zero value of (To, visible quan- 
titative discrepancies are seen when the URMGC and 
HNC/MSApj\/ potential curves are contrasted, or when 
the same is done with the MGC and HNC/MSAppM 
ones. Remarkably, the maximum coincidence between 
the pair of URMGC and HNC/MSApm curves of tpoi'^o) 
(and of i/'/ijp((To)) occurs precisely at the point of zero 
charge. In a similar way, the MGC and HNC/MSAppA/ 
values of ipg (and of ipiHp) converge at ctq = 0, as 
it is known from the studies of the charge-symmetric 
RPM-EDL llOl Il2. ITa. Il6l (for the charge-asymmetric 
case see [Sll l67l iTll l86l l93|l On the other hand, for a 
given value of (To, the differences between the results of 
a Poisson-Boltzmannian theory and the corresponding 
non-punctual HNC/MSA equation for t/jq ^md tJjihp 
grow when size asymmetry is taken into account, e.g. 
at (To = 0.1 C/m2, iijo)^"^^ - (0o)™/*^^^«^*' = 
10.04 mV whereas (V-o)^^^'^^^ - (V^o)^^*^/*^^^^^^ = 
24.11 mV (see Fig. EJa)), and {ipiHpY''^^ - 
{^iHp)"'^^/"^'^''''''' = 10.04 mV whereas 
i.i^iHpy^"''' - (V^mp)™/*^^^-" = 24.12 mV 
(see Fig. EJb)). Also, from Fig. [Sjit is confirmed that 
the potential curves of URMGC approach asymptotically 
to those of MGC when ctq increases, as first pointed out 
by Valleau and Torrie 60]. Contrastingly, the ipo{ao) 
and ipiHp{o'o) curves obtained from HNC/MSAp7\/ 
and HNC/MSAppM exhibit a clear separation for all 
surface charges, which means that counterions do not 
dominate in the SEDL or, else, that the size of the coions 
matters even for large ag. From these observations it 
is infered that: i) an essentially punctual theory (MGC 
or URMGC) is not vahd to describe the SEDL under 
high coupling conditions in which the size of the ions is 
relevant for all their interactions (ion-ion or colloid- ion) , 
except maybe in the neighborhood of ao — 0, and, two 
important points, ii) that including size- asymmetry in 
the model of an EDL is determinant since it exacerbates 
the steric interionic effects previously found in the 
RPM-EDL and iii) that the properties of the EDL are 
not totally determined by the counterions. 

Expectedly, it is found that all the referred phenomena 
detected in monovalent electrolytes are more pronounced 
for divalent ions, as it is evidenced in Figs. IHI^a) andEJb) 
where the associated HNC/MSApm, HNC/MSA rpm, 
URMGC and MGC results for V'o(cto) and ipiMpicro) are 
reported for a 2:2, 0.5 M electrolyte with i?_ = 4.25 A 
and = 8.5 A, around a macrosphere of diameter D = 
160 A and positive charge. For instance, now the discrep- 
ancies between the MEPs of PB and HNC/MSA have 
increased in such a way that even the potentials of zero 
charge, ipo{cro = 0) and i/'/Hp(fo = 0), for URMGC and 
HNC/MSApA/ do not coincide. Once more, URMGC 
has MGC as a limit for ctq — > oo whereas HNC/MSApM 
and HNC/MSAppM do not merge, notwithstanding, this 
time the following very interesting feature shows up: the 
4'iHp{o'o) relationship for HNC/MSAppM exhibits the 
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usual behavior documented in prior studies of positively 
charged colloids in 2:2 solutions 0, l^^l, i-e. it begins 
positive, reaches a maximum and becomes negative only 
at very high surface charges (ctq — 0.344 C/m^), whilst, 
in contrast, ipiHp{<^o) for HNC/MSApA/ behaves differ- 
ently, that is, starts negative at low ctq, then regains "nor- 
mal" positive values, experiences a maximum and, finally, 
becomes negative again at (Tq — 0.22 C/m^, much before 
HNC/MSAijpM does. This new fact could be consequen- 
tial in electrokinetics since the signs of the MEP at the 
IHP and of the electrophoretic mobility of a colloid, /i, 
frequently coincide (due to the identification C = 'ipiHp), 
from which it should be possible to observe experimentally 
a macroion in a multivalent medium having a reentrant 
mobility, i.e. with inversion of fj, at low ctq, a normal 
sign for intermediate charges, and a posterior return to 
a inverted fi for high co- 

A more detailed examination of the size asymme- 
try consequences, already discussed at the level of the 
potential-charge relationship, can be made in terms of the 
SEDL structure. Thence, in the next figures, Figs. [7|and 
|H1 we proceed to analyze the behavior of the ionic pro- 
files of univalent solutions described by the MGC, UR- 
MGC, HNC/MSAflpj,/ and HNC/MSApa/ formahsms. 
To establish a connection with our previous MEP curves, 
the RDFs were obtained for a 1:1, 1 M electrolyte, with 
i?_ = 4.25 A and i?+ = 8.5 A, dissolving a macroion of 
size D = 160 A and surface charges ctq = 0.08356 C/m^ 
(in Fig. 13) and ctq = 0.3004 C/m^ (in Fig. ij). For each 
colloidal charge, the HNC/MSAppM and HNC/MSAp7\/ 
data are compared in the main panels of the figures, 
meanwhile, in the insets the MGC and URMGC plots 
are collated, to evaluate separately the effect of vary- 
ing (To on the distributions gi{r) when size asymmetry 
is incorporated into HNC/MSA or PB theories. From 
the insets it is manifest the almost complete convergence 
between the monotonia URMGC and MGC RDFs for 
growing ag, as it could be awaited from the preced- 
ing results for tpoi^'o) and ipiHp{o'o) H^l- In any case, 
some distinctions between the URMGC and MGC pro- 
files must naturally persist for all ctq (with the only ex- 
ception of the unphysical value cto = oo), due to the 
severe condition g2{r) — for r < (Z? + i?2)/2. Such 
restriction demands a finite discontinuity in g2 ( ) 
for URMGC and thus, for < ctq < oo, impedes a 
perfect coincidence between the RDFs of URMGC and 
MGC in the vicinity of the colloid surface. In other 
words, the total equivalence of URMGC and MGC is 
possible just in the limit of infinite surface charge, since 
only then the infinite and predominant electrostatic re- 
pulsion from the macroparticle will force the colons in 
the MGC description to be completely absent from the 
region {D + Ri)/2 < r < {D + i?2)/2, even if their punc- 
tual nature would allow them to be there. To be more 
explicit, if we define the abbreviations C_ = {D + R-)/2 
and C+ = /2, an inspection of the 1:1 counterion 

contact values, (/_(C_), in TablellTlconfirms the existence 
of small but non-zero differences between the URMGC 



and MGC RDFs close to the macroparticle. Moreover, 
an important piece of information arises from Table UTI 
and the insets of Figs. [7|and|Slfor systems with (Tq > 
and R- < i?+ when tJo is enlarged, namely that the UR- 
MGC counterion contact RDF tends uniformly to that 
of MGC from below (5^«*^<=^C(C_) T 5^^^^(C_)), and 
that the URMGC colon RDF at C+ goes uniformly to 
that of MGC from above {g^^^^^^iC+) i g^^^^{C+)). 
In a more global way, it is noticed therein that, V r, 
gURMGC^^-^ t g^'^'^ir) and gURMGC^^^^ I 5f «^(r) as 

the colloidal charge grows, or, equivalently, that the UR- 
MGC and MGC ionic profiles converge when ctq ^ oo 
but they never cross mutually. 

On the other hand, the comportment of the RDFs for 
HNC/MSApji/ and HNC/MSAppj\/ in the main panels of 
Figs. [3 and IHl contrasts with the PB picture at the insets. 
For our 1:1 electrolyte with counterions smaller than 
colons, the ionic distributions of HNC/MSA in the PM 
and RPM already exhibit a nonmonotonic behavior and 
charge inversion (in the case of RPM these features are 
unapparent due to the scale). Besides, from those graphs 
we arrive to one of the main conclusions of this work, that 
size asymmetry enhances significantly the nonmonoton- 
ical characteristics (e.g. oscillations and charge inver- 
sion) of the SEDL. In relation with the conduct of the 
SEDL structure for high ctq, ihe HNC/MSApm ionic den- 
sities do not have the HNC/MSAhpm ones as a limit 
and, in fact, from Tablenjit is seen that the contact val- 
ues of the counterion RDF in HNC/MSApA^ exceed ap- 
preciably those of HNC/MSAppj\/ even for large values 
of ctq. Additionally, for a constant surface charge density, 
the HNC/MSApniRDFs for counterions and coions fluc- 
tuate more strongly and have a steeper slope than the cor- 
responding to HNC/MSAppM , and, consequently, cross- 
ings between the profiles for PM and RPM do occur. 

For divalent suspensions the RDFs for a 2:2, 0.5 M 
electrolyte, with i?_ = 4.25 A and R+ = 8.5 A, and 
a macrosphere oi D = 160 A are displayed in Fig. |51 
for a surface charge (Tq = 0.08356 C/m^ and in Fig. 
[TUl for CTo = 0.3004 C/m^. Again, in the main pan- 
els HNC/MSApM and HNC/MSAapm are plotted and 
in the insets URMGC and MGC. In general, the situa- 
tion described for the RDFs of monovalent electrolytes 
is repeated here: URMGC and MGC are monotonic, 
whereas HNC/MSApa/ and HNC/MSAijpM osciUate, 
and, for large a^, URMGC converge to MGC with no 
intersections between them and HNC/MSA/jpa/ is not 
the limit of HNC/MSApm. Furthermore, Table HH rati- 
fies that the counterion contact values, .g_(C_) are larger 
for HNC/MSApA/ than for HNC/MSAppA/, and Figs. 
and 1 101 evince that the ionic profiles of the former theory 
wave more intensely than the ones of the latter. Reward- 
ingly, this last pair of structural features of HNC/MSA 
accords with the ARHNC information for a 2:2 planar 
EDL obtained by Greberg and Kjellander | i66i] . 

As a clear example of how the microscopic struc- 
ture determines other properties of the EDL, on the 
grounds of the discussion above we will formulate a 
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rationale for the contrasting high-cro behavior of the 
"tpo^ao) and tpmpi'^o) curves obtained from PB and 
HNC/MSA. Aiming for that, let us first condense some 
of our recent findings in a more convenient form. For 
size-asymmetric SEDL systems with counterions smaller 
than colons, if gi{r]ao) represents an ionic profile cal- 
culated at a given colloidal charge, (i) when (Tq is var- 
ied the URMGC RDFs approach from one side to those 
of MGC, and the separation between the URMGC and 
MGC profiles decreases progressively as ctq is enlarged, 
i.e. there are no crossings between the gf^^''^''^ {r; ao) 
and gf'^^^' {r; (Jo) functions, and the MGC RDFs bound 
those of URMGC when ctq ^ oo, and (ii) the RDFs 
at contact for counterions in HNC/MSApM have larger 
values than in HNC/MSAppA/, and the ionic densities 
of HNC/MSApM have a more pronounced slope than 
the ones of HNC/MSAppM do, i.e. gf (r; ctq) 

and gf^'^^^^^^'^^"{r]aQ) present intersections, and, 
for large ctq, the RDFs of HNC/MSApM are not bounded 
nor have HNC/MSAjipm as a limit. It must be noted 
that in all the previous statements the comparisons be- 
tween the theories are performed at fixed ctq- If we re- 
state explicitely the cro-integral, Eq. for a binary 
electrolyte: 




[g-{t) - g+{t)] t^ dt, (35) 



D/2 



then, for uo 0, this surface charge integral is pro- 
portional to the difference between the areas under the 
r^-functions /_(r;CTo) = g-{r;ao)r'^ and /+(r;cro) = 
g+{r;ao)r'^, for counterions and colons, respectively. In 
such terms, if we consider the data in each of the Figs. 
[3 to IIUI irrespective of the theory, every of the four 
pairs of counterion/coion r^-functions, g-{r;ao)r^ and 
g+{r;a-o)r'^, comprises a constant area. This rephrasing 
of the total electroneutrality condition as a geometrical 
constraint for the RDFs can now be employed to clar- 
ify why the URMGC and MGC potential-charge curves 
converge for ctq — > oo, whereas the corresponding to 
HNC/MSApTif and HNC/MSAppA/ do not, or else, why 
the counterions do not dominate in a description of the 
size-asymmetric SEDL that surpasses the classical PB 
equation (for instance, HNC/MSA). 

To examine the fulfillment of the constant area restric- 
tion as (To is enlarged, in Figs. IlllandlT^the r^- functions, 
f-{r;(Jo) = g-{r;cro)r'^ and /+(r;cro) = g+{r; ao) r'^ , re- 
sulting from the four theories in consideration: URMGC, 
MGC, HNC/MSApM and HNC/MSAppM, are plotted 
(in reduced form) for a 2:2, 0.5 M electrolyte of ionic 
sizes R- — 4.25 A and R+ = 8.5 A, and a macroion of 
diameter Z3 = 160 A bearing two different charge den- 
sities, (To — 0.08356 C/m^ (in the main sections of the 
figures) and ctq — 0.3004 C/m^ (in the insets). For the 
URMGC and MGC data of the lower charge, presented 
in the main panel of Fig. 1111 a comparison between the 
r^-functions for colons immediately shows that some area 



below the URMGC colon function is "lost" in the region 
C- < r < C+, owing to the hard-core condition. On the 
other hard, since the RDFs of URMGC and MGC do not 
intersect, and remembering the PB reciprocity property 
for n : n electrolytes in the region r > C+: 



g^^{r)9^^{r)^l, 



(36) 



it is straightforwardly concluded that the ionic pro- 
files 5^^^'^'^(r; (To) and g^^'^'^'^ {r; ao) are necessar- 
ily located in the region enclosed by g^^'~^'~'' {r;ao) and 
gfGC(^.^^) (i e, Vr, 5™GC(^.^^) < 5*'«^(r;ao) 

and 5+^^^'^'^(r; (To) > g+^^{r;ao)) as the unique way 
in which the colon r^- function of URMGC can recover 
area in r > C+ to assure that the difference functions, 
gd{r;ao)r'^ = [g_(r;(To) ~ g+{r;ao)] r^, of URMGC and 
MGC preserve a constant area. Additionally, when (Tq 
grows (see inset of Fig. MGC colon RDF at 

r = C+ goes down and the colon r^-function of URMGC 
has a minor deficit of area in C_ < r < C+, from which 
the separation between the URMGC and MGC RDFs 
diminishes and, eventually, disappears at ctq = oo. By 
contrast, in the HNC/MSA case (Fig. there is also 
a lost of area in the zone C_ < r < C+ for the colon 
r^-function of HNC/MSApM, however the possibilty of 
crossings between the ionic profiles of HNC/MSApM 
and HNC/MSAppM, stemming from the enhancement 
of the non-monotonicity of the RDFs caused by the 
size-asymmetry, does not obligate the coincidence of 
^HNC/MSA,., (r; (To) and gf^c/MSA,.., ^^^^^ 

(To , even if the difference functions of both theories have 
the same area. It is from these reasonings that the con- 
vergence, for (To — * oo, between the structural proper- 
ties of URMGC and MGC and the non-convergence for 
those of HNC/MSApM and HNC/MSAppA/ can be un- 
derstood. A similar analysis done for 1:1 systems (not 
shown) confirms all the arguments exposed. 

Summing up, when the value of ao rises, the dif- 
ference function gd{r;cFo)r^ = [g-{r;ao) - g+{r;aa)] 
for each of the MGC, URMGC, HNC/MSAppM and 
HNC/MSApAf approximations can accomplish the con- 
straint of a constant area in two ways. On the one 
hand, when (Tq oo, the r^-integrals (or second mo- 
ments) of 5y^*^'^'^(?'; (To) and g^^''^'^ {r; (Jq) remain ex- 
actly equal whereas all the rest of integrals of such 
functions (including the ipo- and "i/jj/fp-integrals; cf. 
Eqs. f (|27() and l^'29\i ) become progressively alike, since 
the RDFs of URMGC and MGC approach asymp- 
totically. On the other, when ao increases, the 
pair 5f ^^/^^^^p- (,; and g^^c/MSA...., ^^^^ 

the same r^-integral but other integrals (such as the 
Ipo- and "i/j/^rp-integrals) become different, since the 
HNC/MSApM and HNC/MSAppM ionic profiles do 
not merge. This cogently explains why the ■0o(o'o) 
and tpiHp{o'o) curves for URMGC and MGC con- 
verge for large ao and those for HNC/MSApM and 
HNC/MSAppA/ remain separated. 
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To complete this study of the PM-SEDL we switch 
our attention to the MEP and CRC profiles, ip{r) and 
Q*(r), respectively, to gain more insight into the com- 
portment of the diffuse EDL in terms of its neutraliza- 
tion (or charge screening) capacity. In Figs. 1131 and [T^ 
the MEP and CRC functions corresponding to the 2:2 
data in Figs. El and ^| are portrayed, respectively. We 
have preferred to graph the divalent and not the univa- 
lent results since the effects to be discussed have qual- 
itative similarities in the 2:2 and 1:1 systems but are 
more marked in the former case. In respect to the tp{r) 
profiles, the URMGC and MGC MEPs included in the 
inset of Fig. E| confirm the monotonic character of the 
PB solutions and the convergence between URMGC and 
MGC for (To ^ oo. In turn, all the HNC/MSA MEP 
profiles in the main panel of Fig. 1131 oscillate and present 
ample regions of inverted potentials or even show an al- 
ternation of signs (see, for example, HNC/MSApj\/ for 
ctq = 0.3004 C/m^). However, the most salient fea- 
ture in both the PB and HNC/MSA results of tpir) is 
that, in the proximities of the colloid, an SEDL de- 
scription in which size asymmetry is assumed exhibits 
lower MEP values compared to the corresponding poten- 
tials of the size-symmetric version of the same theory, 
i.e. 'i{j"^^/'^^^''^^^'{r;cro) < ^hnc/msa^pm (^j..^^) and 
^URMGC^^.^^-^ < ^^(^(^{r-ao), for {r - D/2) < 7.5 A. 
A complementary and more quantitative estimation of 
this can be obtained from the values collected in Ta- 
ble mil The observed behavior of the ip{r) functions 
close to the surface implies that a size-asymmetric dif- 
fuse EDL (either of HNC/MSApa/ or URMGC) neu- 
tralizes more efficiently the native macroion charge than 
that of a size-symmetric system (of HNC/MSA^pm 
or MGC). This generalized dominance of the neutral- 
ization power of an unequally-sized electrolyte is re- 
markable since in the PB approaches there are less 
counterions close to the macroparticle in URMGC 
than in MGC, in other words, the counterion contact 
RDFs comply with g^^^^^^ {C-;ao) < gi^G^(C_; ctq), 
whilst in the HNC/MSA theories the contact number 
of counterions in HNC/MSAp7\/ exceeds the predic- 
tion of HNC/MSAppM, i.e. .g!^^^/^^^^^"(C_;(To) > 
g^^(^/^tSARPM ^(j .^^^y Such "anomaly" can be jus- 
tified nicely by the obligated absence of coions in the 
zone C- <r <C+ for the HNC/MSApAf and URMGC 
theories, which, irrespective of the PB or HNC/MSA 
approach, allows the counterions to compensate the 
colloidal charge in a better way than it is done in 
HNC/MSAppM or MGC, theories where the coions can 
be indeed present in C_ < r < C+ due to the condition 
— This, however, does not explain the addi- 

tional neutralization power of HNC/MSApM over UR- 
MGC. For HNC/MSApM the reason is the increase of 
the excluded volume, which in turns promotes the coions 
and counterions adsorption, due to entropic effects, as 
it has been discussed by Jimenez- Angeles and Lozada- 
Cassou . These last effects also explicate the increase 
of the counterion adsorption of the HNC/MSApj\/ over 



HNC/MSAppA/ and, consequently, the enhancement of 
the charge screening in the PM-SEDL. 

The behavior of the neutralization capacity detected 
via the mean electrostatic potential profiles fits very well 
with the information provided by the cumulative reduced 
charges. For instance, in the inset of Fig. E| it can 
be verified that the Q*{r) functions for URMGC and 
MGC are monotonic and do not display any change 
of sign, contrasting with the main panel of the same 
figure where HNC/MSApa/ and HNC/MSAppA/ have 
CRC curves with oscillations and appreciable charge re- 
version. Additionally, a genuine and clearer manifes- 
tation of the enhanced neutralization capacity of size- 
asymmetric models over the symmetrical ones can be 
gathered directly from the full set of CRC curves in Fig. 
1141 There, it is visible the faster decay in the cumu- 
lative charges corresponding to HNC/MSApa/ and UR- 
MGC with respect to those of HNC/MSAppM and MGC, 

i.e. Q*"^^'^'SAPM^^.„^) < Q*HNC/MSAnPM^^.„^) 

and Q*^^*^'5'^(r;cro) < Q*^'^^ {r^a^), for (r - D/2) < 
10 A. 

Notably, for a given cto, among all the reviewed the- 
ories, the HNC/MSApAf formalism is always associated 
to the strongest fluctuations and to the more accentu- 
ated sign inversions in ip{r) and to the largest neutral- 
ization and charge reversal effects in Q*{r). This fact, 
together with the rest of findings revealed by the anal- 
ysis of the ip{r) and Q*{r) profiles, conforms an extra 
and compelling testimony to the relevance of using size- 
asymmetric models when a more faithful description of 
colloidal suspensions is sought. 

To end, we want to point out that, for our n : n sys- 
tems, at the distances where the RDFs show a charge 
inversion the CRCs go through a minimum or maximum, 
as it was noted previously by Wang et al. in the RPM- 
EDL for cylindrical geometry This signifies that, for 
a valence-symmetric RPM or PM electrolyte, charge in- 
version implies charge reversal. What is more, from basic 
electrostatics for the SEDL, it is possible to prove that 
^^^7^ = ~^Tr^ (similar equations for other geometries 
are given in Refs. [I^jEIDj thus, the distances at which 
the CRC is zero correspond precisely to those where the 
mean electrostatic potential presents extremum values, 
whence, the existence of oscillations in Q*{r) should not 
be enough to produce fluctuations in V'('')j i-S. a change 
of sign in the CRC is indeed required. 



IV. CONCLUSIONS 

In this work a survey of the size-asymmetric spherical 
electrical double layer in the primitive model was per- 
formed by using the HNC/MSA integral equation. Af- 
ter the correctness of the HNC/MSA description was at- 
tested by comparing our numerical results with simula- 
tions and previous ARHNC data (in the planar limit), 
we have carried out a study of the structural and charge- 
potential relationship for archetypal cases of 1:1 and 2:2 
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systems with a size-asymmetry ratio of 2. To assess 
the importance of the size-asymmetry effects, we col- 
lated our predictions with those of HNC/MSA in the 
RPM, and also with the corresponding to the classical 
Poisson-Boltzmann approaches: URMGC and MGC, to 
exhibit the notable quantitative and qualitative discrep- 
ances existing between punctual and non-punctual EDL 
theories. The main conclusion of this paper is that size- 
asymmetry is an essential improvement to consider in 
studies of the SEDL since heightens ionic size correla- 
tion effects (e.g. nonmonotonicity, charge inversion and 
charge reversal) already seen in the RPM case and, conse- 
quently, unveils new and interesting phenomenology ab- 
sent in size-symmetric systems. As a conspicuous exam- 
ple of these novel effects, it is evidenced that counterions 
do not always dominate in the PM-SEDL, as it is sug- 
gested by the PB point of view, and the wherefores of 
this fact are traced back to the detailed characteristics of 
the EDL structure. In particular, the contrasting high-cro 
behavior displayed by the PB and HNC /MSA potential- 
charge relationships, i.e. the convergence, for ctq oo, 
between the URMGC and MGC VoCo-q) and iIjihp{<Jo) 



curves and the non- merging between the HNC/MSApM 
and HNC/MSAtjpa/ ones, was elucidated by the follow- 
ing rational: for a system with fixed parameters, when 
CTo increases the RDFs of URMGC and MGC, and all 
their integrals and associated EDL properties (e.g. ipo 
and ipiHp), go similar, whereas for HNC/MSApM and 
HNC/MSAjiPM , even if their RDFs can have the same 
second moment, in general, their corresponding shapes 
and ensuing properties (ipo and/or i/jihp) could be com- 
pletely different. 

In forthcoming publications some other appealing phe- 
nomena occurring in PM-SEDL, as overcharging lip and 
the advent of anomalous differential capacities j93|, will 
be addressed by means of our integral equation approach. 
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APPENDIX A: EXPLICIT MSA BULK DIRECT 
CORRELATION FUNCTIONS 

The direct correlation functions of a PM bulk elec- 
trolyte in the mean spherical approximation have been 
obtained analytically by Blum and Hiroike 

for the general case of n species with arbitrary size and 
charge. The only restriction is the electroneutrality con- 
dition: 



1~ «(/0li??+P2i? 




and 



(A9) 



eksT 



(AlO) 



^ PiZ, = 0. 



with pi and p2 being the numerical densities of the two 
(Al) species. 

If we define 



The MSA direct correlation functions can be expressed 



i?2 — ^1 



(All) 



iis) 



(A2) 



where c^j{s) is the hard-sphere contribution and c^j"^^ 
is the electrostatic part. 



and 



Ri2 — 



Ri + i?2 



(A12) 



1. Electrostatic contribution 

Following the Baxter method Blum found that the 
excess properties of an electrolyte can be written in terms 
of the parameter T. When a binary mixture of hard 
charged spheres of diameters Ri and R2 and valences zi 
and Z2 is considered, this parameter can be obtained solv- 
ing the trascendental algebraic equation (it is assumed for 
definiteness that R2 > Ri). 



_ ZjZjC 



(A13) 



then cf^'^{s) ~ 1,2) is given exphcitly by 



Ci2''(s) = /3o,0<s<A, 



(AM) 



= t:loD{T), 



(A3) cf2^''{s)^aas-'^+ai + a2S + a3S^,X<s<Ri2, (A15) 



where 



D{T) = piXi\T) + P2X2^{T), (A4) 



4-(5) = -Pi2S-\s>i?i2, (A16) 



where 



^i(r) = 



^1 , RlxiT) 



^2(r) 



1 + ri?i 1 + TRi ' 



22 , i?ix(r) 



1 ri?2 1 + TR2 ' 



(A5) 



(A6) 



3(n-ri?i) 



Zl(i?2--Rl)\ f Z1Z2T 



T 



i + ri?2 



, (A17) 



~c {piRiZijl + ri?i)-i + P2R2Z2{1 + TR2)-^} ao 
1 + c{pii?3(l + ri?i)-i + p2i?i(l + TR2)-^} ' 

(A7) 



T = (i + ri?i)(i + ri?2), (A8) 



uj{R2-Ri)^ 
16T 



4:{Rl + Rl) - AT'^RIrI 



-{Ri-R2rT 



4(zi + Z2) 



X + 4ziZ2r' 



(A18) 
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where 



Zl(R2 - Rl) + Z2(Rl - R2) 

-Rl 



[Rl + R2)ziZ2T' 



X - 2TziZ2 

(i + ri?2)(ri?i -2) 



Rl 



(i + ri?i)(ri?2 -2) 



X 



(A19) 



Jo = 



{I + TR2Y 



- 2Vzl - 2R,2zIt'^ 



';Rl{TR2-2){l + TR2) 



(A29) 



^2 = Y 



2^1 + ^2 X + r Z1Z2 

nl^ I T^'2 r?2 7;)2 



2 , n2 T-2D2n2 (i?l-i?2)^T 



and 



as 



wx 



cf^'^is) = io + «is + i2s\ < s < (A22) 



s) = -Piis \ S > i?l 



(A23) 



with 



a; 



2t^2 



(i + ri?i)2 



2rzf - 2i?izfr 



^i??(ri?i-2)(i + ri?i^ 



(A24) 



^l 



(l + TRi) 
^ [2zi] X + T^zf 



— (^[2Rj^irRir] 



and 



«2 



(A25) 



(A26) 



Finally, 



C22" is) = Jo + Ji s + J2s' , < s < i?2 , (A27) 



cg'%s) = -P22S-\s> R2 



(A28) 



Ji = 



(i + ri?2)2 

[2Rl~{TRir]x 



[2z2] X + r'zf + 



(A30) 



(A20) and 



(A21) 



J2 



WX 



(A31) 



2. Hard sphere contribution 



The exact solution for the Percus-Yevick equation of 
a hard spheres mixture with n species was analyticaly 
obtained by Lebowitz jo^. Following the notation of 
Lebowitz, we define for a binary mixture 



e = - {piRl + P2RI) 



(A32) 



ai = (1 - <l+^ + e+ qRKpi + P2)(l + 20 



-{R2-Rlfp2 



-R2PiRl + {Ri+R2) 
6 



:RiR2{piR\ + P2RI) 



[pi + P2)[l + e + - ^PlP2iR2 - Rlf 



{Rl + R2) + -RiR2{piRl + P2RI] 



(A33) 



02 = (1 - ir' ii+^+e + -^RiiPi + P2)(i + 20 



-{R2-Rifpi 



7-R1/92-R2 + (^1 + -^2 



:RiR2{piR\ + P2Rl) 
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(pi + P2)[l + ^ + e]- ^PlP2iR2 - Rlf 



{Ri + R2) + -RiR2{piRl + P2RI) 

D 



(A34) 



where 



7o = A2(6-3(iA"), 



(A47) 



= [Pia\ + Pia2] , 



(A35) 



gu{Ri 



2 4' 



l + | + ^P2i?^(^l-^i2)}(l-0-', 



(A36) 



2 4' 



-i^l)|(l-0-^ 

(A37) 



71 = 8A-^rf - 2A6 



and 



72 = 6 - SA^d. 
Finally, if we define the constants 



(A48) 



(A49) 



Ai=io-ai; A2=ii-bi; A3 = 12 - d; (A50) 



9i2{Ri2) = [R2gii{Ri) + Ri922{R2)]{Ri+R2) \ Bi=jo-a2; ^2=^1- 62; B3=j2-d; (A51) 

(A38) 



&i = -4 [4pi-R?5ii^(^ii) +P2(i?l +i?2)V2^(iii2)] , 

(A39) 



Ei=Po-ai; Fi=ao-7o; (A52) 



62 = -J [^P2Ri922^{R2) + Pl{Rl + R2f9l2^{Rl2)] 



and 



(A40) 



b = -- 



piRigii{Ri) + 



P2R2922{R2) 



{Ri+R2)gi2{Ri2). (A41) 



The c'^f{s) = 1,2) are then given by 



-F2 = "1 -fli -71; -^3 = 0:2-72; ^4=0:3-6?; (A53) 
we can write the bulk correlation functions as 



Cii(s) = c^(s) + d^'is) = c!-\s) + c^fis), (A54) 



where 



c\f{s) = Piis-^+Ai + A2S + A3s^,0< s < Ri, (A55) 



cf/(s) = -ai - biS - ds^ s<Ri,i=l, 2, (A42) 



cif (s)=0,s>i?i. 



(A56) 



c':f{s)=0,s>Ri,i = l,2, (A43) c^^'{s) = P22S-' + B, + B2S + B3s\0 < s < R2, (A57) 



cfl(s) = 4t{s) = -ai, < s < A, (A44) 



4^*(s)=0,s>i?2, 



(A58) 



c^2(s)=4^(s) = -ai-7os-' 
— 71 — 72S — ds^, A < S < iii2. 



(A45) 



c™*(s) = 4f (s) = Pi2S-^ + i^i, < s < A, (A59) 



4^is)=4l{s)=0,s> R,2, (A46) 



mi / \ int / \ j-t — 1 , 

-FlS"^ +F2 +i^^3S + ^4S^A < S < i?i2, (A60) 
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where 



cif(s) = <(s) = 0,s>i?i2, (A61) 



4f{s)^-Pns-\s>0, 



4f{s) = -P22S-\s>0 



(A62) 



(A63) and 



Kha{r, t) = 27rgi {1)1^ [PnJiRi, -l,r,t) + 
AiJ{Ri,0,r, t) + A2J{Rul,r, t) + 
A3J(i?i,3,r,t) ) (B7) 



and 



clf{s)=clf{s) = ^P,2S-\s>Q. 



(A64) 



Khh{r,t) ^ -2^Piigi{t)-[r + \ r - t\ ). (B8) 



If 



APPENDIX B: EXPLICIT FORM OF THE h{r) 

AND Hi{r) TERMS IN THE HNC/MSA 
INTEGRAL EQUATIONS FOR THE PM-SEDL 



D Ri 



< r < oo, 



(B9) 



Due to the radial symmetry of the system, the angular 
dependence {9,(j)) of the equations ((TTJ-I^^l), l|T^ - l|T^ 
and (|21|) - (|24|l can be resolved since the angular integrals 
can be performed analytically. 

The final form of the Ii{r) and Hi{r) terms can be 
given if we define 



then 



hir) = 47r 



PiiRl AiRl A2Rt AsRf 
2 ^ 3 ^ 4 ^ 6 



If 



(BIO) 



J{xfJ,r,t) = 



1 



(l + 2)rt 



4+2- |r-t|'+2 



and 



(Bl) 



then 



D Ri D R2 , , , 

— + — <r<— + — + A, (Bll 
2 2--2 2 ' ^ ' 



J'{xi,Xf,l,r, t) 



[l + 2)rt 



J+2 _ 1+2 
Xf Xi 



Thus, if 



we have 



D Ri D 3^ 

h— <r< h -Ri, 

2 2 - - 2 2 ^' 



(B2) 



(B3) 



h{r) = 



r+\ 



2^2 



Kl3a{r, t)dt 



r+Ri 



r+X 



Khb{r,t)dt 



D I «2 
2 2 



Khc{r, t)dt, 



(B12) 



if 



hir) 
and, if 



r+Ri 



KIia{r,t)dt+ / KIib{r,t)dt, 



f + f + A<.<f + f + (B13) 



(B4) then 



then 



hir) 



D 3„ 

— + -Ri <r <(X), 



(B5) 



pr+Ri poo 

/ Kha{r,t)dt+ KIib{r,t)dt, (B6) 



hir) 



r — X pr-\-X 

Khbir,t)dt+ / Khair,t)dt 



2^2 

+ / Khbir,t)dt 

Jr+X 



r-X 



D I -"2 
2^2 



Kl3cir,t)dt, 



(B14) 



and, if 
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then 



^ ^2 „ 

+ ^ + Ri2 < r < oo, 



(B15) 



then 



pr — X fr-\-X 

h{r)= KIsb{r,t)dt+ KIsa{r,t)dt 

Jr—R-in Jr — \ 



h — Ri2 Jr — \ 

fr+Ri2 roo 

+ I Khb{r,t)dt+ / Kl3c{r,t)dt, 



r+X 



D . «2 
2 2 



(B16) 



where 



Khair, t) = 2ng2it)t'^ (^Pi2 J(A, -1, r, t) + 

E^J{X,0,r,t) + (Pi2 + Fi)J'{X,Ri2,-l,r,t) + 
Fa J'(A, Ri2, 0, r, t) + FzJ'{X, R12, 1, r, t) + 



F4J'(A,i?i2,3,r,i) 



(B17) 



Khb{r, t) = 2ng2{t)t'' (^Fa J(iii2, 0, r, t) + 
F3 J(i?i2, 1, r, t) + F4 J(i?i2, 3, r, + 
(Pi2 + i^i)J(i?i2,-l,r-,t)') (B18) 



and 



Khc{r, t) = -27rPi2ff2(i)- [r + t-\r-t \ ) . (B19) 



If 



D Ri 
2+T 



< r < 00, 



(B20) 



then 



Hi 



(r)= / ^ KHib{r,t)dt+ / KHia{r,t)dt 



2^2 



r-A 



/■r+i?i2 /'oo 

+ / KHib{r,t)dt+ / KHic{r,t)dt, 
Jr+X J§ + ^ 



and, if 



Y + + -K12 < r < 00, 



(B23) 



(B24) 



then 



nr — X pr-\-X 

Hi{r)= / KHib{r,t)dt+ / KHia{r,t)dt 

Jr—Ri2 J r—X 

rr-\-Ri2 roo 

+ / KHib(r,t)dt+ / KHic{r,t)dt, 

Jr+X + ^ 

(B25) 



where 



/<:ifia(r, i) = 27751 (^Pi2 J(A, -1, r, t) + 

EiJ{X, 0, r, t) + (P12 + Fi) J'(A, J?i2, -1, r, i) + 
F2 J'(A, iii2, 0, r, + F3 J'(A, iii2, 1, r, + 



F4J'(A,Pi2,3,r,t)^, 



(B26) 



KHib{r,t) = 2T,gr{t)e\F2J{Ri2,(i,r,t) + 
F3 J(i?i2, l,r,t) + F4 J(i?i2, 3, r, + 
(Fi2+Fi)J(i?i2,-l,r,i) ) (B27) 



and 



/ \2 \3 

/4(r) =47rfPi2y +Fiy + 



(Fi + Pr. 



2 

+ F4- 



+ F2- 



+ 



A<^) 



(B21) 



If 



KHic{r, t) 



If 



-27rPi25i(i)- r + i- | r-t 
r 



D R2 

_ + _<r<oc, 



then 



(B28) 



(B29) 



H2{r)=h{r). 



(B30) 
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If BiJ{R2,0,r,t) + B2J{R2,l,r,t) + 

BsJiR2,3,r,t)) (B35) 

then and 

pr+R2 poo 

H^{r)= KH3a{r,t)dt+ KHzb{r,t)dt, t( \ 

(B32) ^ V / 



and, if 



Finally, if 



^ + ^i?2 < r < 00, (B33) 



then ^ + ?^<r<oo, (B37) 



H3{r)= KH3a{r,t)dt+ KH3b{r,t)dt, then 

Jr-R2 J^ + ^ 

(B34) 



where 



6 

KH3a{r, t) = 2tt g2{t)e I F22^(i?2, -1, r, i) + (B38) 



Ri/2 




FIG. 1: Schematic representation of the model. 




r'(A) 

FIG. 2: SEDL radial distribution functions for a 2:2, 0.5 M 
electrolyte around a colloid of diameter and surface charge 
density _D = 10 A and ao = 0.407 C/m^, respectively. The 
ionic species have diameters R- = 4.5 A and R+ = 8.5 A. 
In the main panel the circles and the continuous lines cor- 
respond to Monte Carlo data and HNC/MSApAf results, re- 
spectively. In the inset the same system is considered, with 
the dashed and continuous lines corresponding to molecular 
dynamics simulations and to the HNC/MSApAf theory, re- 
spectively. The distance r' is measured from the macroparti- 
cle's surface. 



20 




r'(A) 



FIG. 3: SEDL radial distribution functions for a 2:2, 0.5 M 
electrolyte around a colloid of diameter and surface charge 
density D = 10 A and ctq = 0.407 C/m^, respectively. The 
ionic species have diameters R- = 8.5 A and R+ = 4.25 A. 
The symbols and curves have the same meaning as in Fig. |21 




r'(A) 



FIG. 4: HNC/MSAPA/ and ARHNC radial distribution func- 
tions for monovalent and divalent 1 M electrolytes. For 
HNC/MSApA/ (continuous lines) the ionic distributions are 
around a very large colloid of diameter D = 1000 x R- and 
surface charge density ctq = 0.267 C/m^. For the ARHNC 
theory (with symbols) the profiles are next to a wall of the 
same charge ao = 0.267 C/ra^. In the main panel the 1:1 
ionic species have diameters R- = 4.25 A and R+ = 8.5 A 
and in the inset the 2:2 ionic diameters are R- = 4.25 A and 
i?+ = 6.375 A. The distance r' is measured from the charged 
surface. 




a„(C/m) 



FIG. 5: SEDL mean electrostatic potential at; [a] the surface 
of the colloid and [b] the inner Helmholtz plane, as a function 

of the surface charge density uq, for a 1:1, 1 M electrolyte 
with R- = 4.25 A and R+ — 8.5 A around a macroion of 
diameter D = 160 A and non-negative surface charge den- 
sity. The solid and dotted lines stand for HNC/MSApm 
and HNC/MSAiiPM, and the dashed and dot-dashed lines 
correspond to URMGC and MGC, respectively. The ion- 
colloid closest approach distances arc: d- — {D + R-)/2 and 
d+ = {D + R+)/2 for HNC/MSApm and URMGC, and d- = 
d+ = {D + R-)/2 for HNC/MSAppm and MGC. The ion-ion 
closest approach distances are: d-+ = d+- = {R- + R+)/2 
for HNC/MSApM, d-+ = d+- = R- for HNC/MSAhpm, 
and d_ + = d+ - = for URMGC and MGC. 




TABLE I: Parameters of the 2:2, 0.5 electrolytes used in 
the Monte Carlo and molecular dynamics simulations of the 
asymmetric SEDL. A macroion with diameter D = 10 A and 
surface charge density ao = 0.407 C/m^ is employed in all the 
runs. 



Run Method zi Z2 Ri (A) R2 (A) iVi N2 L 



A 


MC 


-2 


+2 


4.25 


8.50 


1308 


1304 


38.35 Ri 


B 


DM 


-2 


+2 


4.25 


8.50 


2894 


2890 


50.00 Ri 


C 


MC 


+2 


-2 


4.25 


8.50 


1304 


1308 


38.35 Ri 


D 


DM 


+2 


-2 


4.25 


8.50 


2890 


2894 


50.00 Ri 
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FIG. 7: SEDL radial distribution functions for a 1:1, 1 M 
electrolyte with R- = 4.25 A and R+ = 8.5 A next to a 
colloid of diameter D = 160 A and surface charge density ag — 
0.08356 C/m^ In the main panel the solid and dashed lines 
stand for HNC/MSAhpm and HNC/MSApm, respectively, 
whereas in the inset the solid and dashed lines correspond to 
MGC and URMGC, respectively. The colloid-ion and ion-ion 
closest approach distances used in each theory are the same 
as in Fig. |3 Here (and in the rest of figures) the distance r' 
is measured from the macroparticle's surface. 



TABLE 11: SEDL values of the radial distribution functions at C- = {D + R-)/2 and C+ = {D + R+)/2 for 1:1, 1 M and 2:2, 
0.5 M size-asymmetric electrolytes around a colloid of diameter D = 160 A and variable ctq, obtained from the HNC/MSApM, 
HNC/MSAppA/, URMGC and MGC theories. The ionic species have diameters R- — 4.25 A for the counterions and R+ — 8.5 
A for the colons. The colloid-ion and ion-ion closest approach distances used in each theory are the same as in Fig. |3 The 
surface charge density ctq is in C/m^. 



Electrolyte 


o-Q 


^HNC/MSApj,j 


^HNC/MSAj,jv, 


^HNC/MSAHf>A/ 


HNC/MSA^Pj,,f 

9+ 


gURMGC 


URMGC 
9+ 


gMGC 


gMGG 


1:1 


0.08356 


5.87 


1.6 


3.7 


0.3 


3.1 


0.6 


3.4 


0.3 


1:1 


0.3004 


28.8 


0.06 


25.4 


0.01 


24.4 


0.3 


24.6 


0.04 


2:2 


0.08356 


6.2 


0.5 


5.6 


0.2 


4.9 


0.6 


5.3 


0.2 


2:2 


0.3004 


50.7 


0.03 


49.7 


0.009 


47.8 


0.3 


48.1 


0.02 
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FIG. 8: The same as in Fig. |3but for a surface charge density 
ao = 0.3004 C/m^ 



TABLE III: SEDL mean electrostatic potentials at the macroion's surface (tpo) and at the inner Helmholtz plane {ijiiHp) for 
the HNC/MSApM, HNC/MSAapm, URMGC and MGC theories. The parameters are the same as in Table lUl The surface 
charge density (tq is in C /m? and the MEPs in mV . 



Electrolyte 




^HNC/MSAj,jv, 


,HNC/MSApjv^ 
VlHP 


^HNC/MSAhpm 


, HNC/MSAjjpju 
VlHP 


j,URMGC 

Vo 


j,URMGC 
WlHP 




J.MGC 
WlHP 


1:1 


0.08356 


33.4 


8.5 


47.8 


22.9 


54.1 


29.2 


56.4 


31.5 


1:1 


0.3004 


124.3 


34.8 


147.4 


58.0 


171.5 


82.0 


171.7 


82.2 


2:2 


0.08356 


26.7 


1.8 


36.2 


11.3 


45.3 


20.4 


46.3 


21.4 


2:2 


0.3004 


82.6 


-6.9 


94.2 


4.8 


139.2 


49.7 


139.2 


49.7 
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FIG. 9: SEDL radial distribution functions for a 2:2, 0.5 M 
electrolyte with R- = 4.25 A and R+ = 8.5 A next to a 
colloid of diameter D = 160 A and surface charge density cro = 
0.08356 C/m^ In the main panel the solid and dashed lines 
stand for HNC/MSAhpa/ and HNC/MSApm, respectively, 
whereas in the inset the solid and dashed lines correspond to 
MGC and URMGC, respectively. The colloid-ion and ion-ion 
closest approach distances used in each theory are the same 
as in Fig. El 




FIG. 10: The same as in Fig. |5| but for a surface charge 
density ctq = 0.3004 C/m^ 
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r' (A) 

FIG. 11: SEDL reduced profiles /*(r') = }{{D/2) +r')/Cl 
for a 2:2, 0.5 M with _R_ = 4.25 A and R+ = 8.5 A around 
a colloid of diameter D — 160 A . The solid and dashed lines 
are for MGC and URMGC, respectively. In the main panel 
the surface charge density is ctq = 0.08356 C/m^, and in the 
inset is ao = 0.3004 C/m^ The colloid -ion and ion-ion closest 
approach distances used in each theory are the same as in Fig. 
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r' (A) 



FIG. 12: SEDL reduced profiles /*(r') = f{{D/2) +r')/Cl 
for a 2:2, 0.5 M with _R_ = 4.25 A and R+ = 8.5 A around a 
colloid of diameter D = 160 A. The solid and dashed lines are 
for HNC/MSAhpm and HNC/MSApa/, respectively. In the 
main panel the surface charge density is ao = 0.08356 C/m^, 
and in the inset is (tq = 0.3004 C/m'^. The colloid-ion and 
ion-ion closest approach distances used in each theory are the 
same as in Fig. |3 
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FIG. 13: SEDL mean electrostatic potential as a function 
of the distance to the macroparticle's surface for a 2:2, 0.5 
M with R- = 4.25 A and R+ = 8.5 A around a colloid of 
diameter D — 160 A . In the main panel the dotted line is for 
HNC/MSApAf and ao = 0.08356 C/m^ the solid line is for 
HNC/MSApM and ao = 0.3004 C/m^ the dot-dashed is for 
HNC/MSAflPM and ao = 0.08356 C/m^ and the dashed line 
is for HNC/MSAflPM and ao = 0.3004 C/m^ In the inset the 
dotted line is for URMGC and ao = 0.08356 C/m^, the solid 
line is for URMGC and ao = 0.3004 C/m^, the dot-dashed 
is for MGC and ao — 0.08356 C/m^, and the dashed line is 
for MGC and ao = 0.3004 C/m^ The colloid-ion and ion- ion 
closest approach distances used in each theory are the same 
as in Fig. |^ 
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FIG. 14: SEDL cumulative reduced charge as a function of 
the distance to the macroparticle's surface for a 2:2, 0.5 M 
with R- = 4.25 A and R+ = 8.5 A around a colloid of di- 
ameter D = 160 A . In the main panel the dotted line is for 
HNC/MSApM and ctq = 0.08356 C/m^ the solid line is for 
HNC/MSApM and ao = 0.3004 C/m^, the dot-dashed is for 
HNC/MSAflPA/ and ao = 0.08356 C/m^ and the dashed line 
is for HNC/MSAhpa/ and ao = 0.3004 C/m^. In the inset the 
dotted line is for URMGC and ao = 0.08356 C/m^ the solid 
line is for URMGC and ao = 0.3004 C/m^, the dot-dashed 
is for MGC and ao = 0.08356 C/m^, and the dashed line is 
for MGC and ao = 0.3004 C/m^ The colloid -ion and ion-ion 
closest approach distances used in each theory are the same 
as in Fig. |S1 



